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Abstract 

We incorporate the effect of lattice thermal vibrations into the Glauber theory description 
of particle- and nucleus- crystal Coulomb interactions at high-energy The allowance for the 
lattice thermal vibrations is shown to produce strong absorption effect: the phase shift 
function of the multiple-diffraction scattering on a chain of N identical atoms acquires large 
imaginary part and the radius of the absorption region in the impact parameter plane grows 
logarithmically with N. Consequences of this observation for the elastic and quasi-elastic 
Coulomb scattering are discussed. Practically interesting example of the coherent Coulomb 
excitation of ultra-relativistic particles and nuclei passing through the crystal is considered 
in detail. 
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1 Introduction 



In this paper we develop the Glauber theory [1] description of the absorption phenomenon 
in coherent particle- and nucleus-crystal Coulomb interactions at high-energy. 

As is well known, the multi-loop corrections generate the imaginary part of the scattering 
amplitude even if the tree-level amplitude is purely real. For example, the purely real Born 
amplitude of the high-energy Coulomb scattering in crystal acquires the imaginary part due 
to the multiple scattering (MS) effects. However, in the widely used static/frozen lattice 
approximation (SL approximation) the account of rescatterings alters only the overall real 
phase of the full amplitude thus producing no absorption effect . The latter is related to 
the creation and annihilation of excited intermediate states of crystal and as such manifests 
itself only beyond the SL approximation (for the analysis of elastic scattering based on the 
SL approximation see [2]). Indeed, the amplitude of small-angle elastic scattering on a chain 
of N identical atoms in the impact parameter representation equals 

1 - (S(b)), 

where the scattering matrix placed between the ground states of crystal is 

(S(b)) = (exp[i X (6)]> 

with the purely real phase shift function x(p) = YljLiXjiP)- I n the SL approximation 

(exp[ix(&)]) ~ exp[ix(&)]. 

Therefore, only with the allowance for the lattice thermal vibrations the Coulomb phase 
shift function gets, in general, non-vanishing imaginary part which is interpreted as an 
absorption effect. The imaginary part appears only as the second order perturbation, 

- ^(x 2 ) - m 

But the strength of the effect is oc (3 2 N, where (3 is the coupling constant. For the coherent 
scattering of relativistic nuclei (the electric charge Z\) on the chain of N atoms (the atomic 
number Z 2 ) in a crystal the effective coupling (3 = 2aZ\Z 2 is strong and the absorption 
effect is strong as well. The absorption is strong for impact parameters, b, smaller than 
some characteristic value b a oc log(/3iV), and vanishes toward the region of larger b. This 
phenomenon provides a natural ultra-violet regulator of the theory and enables, in partic- 
ular, consistent calculation of the coherent elastic scattering cross section. The latter is 
calculated and turns out to be equal to the one half of the total cross section. As we shall 
see, the absorption effect is of prime importance also for quantitative understanding of the 
phenomenon of the coherent Coulomb excitation of relativistic particles and nuclei passing 
through the crystal. A consistent description of this phenomenon is the goal of our paper. 

The outline of the paper is as follows. We start with the well known example of the 
coherent Coulomb elastic scattering of charged particle/nucleus by a linear chain of N 
identical atoms in a crystal target (Sec. 2). We derive the scattering matrix with absorption 
and calculate the cross section of the coherent elastic scattering, a e \ (Sec. 3) and the cross 
section of the incoherent excitation and break-up of the target, <jQ e i (Sec. 4). We find that 
in the large- N limit a e i ~ o-Q e i ~ \a to t (Sec. 5). In Sec. 6 we discuss the coherent Coulomb 
excitation of ultra-relativistic particles and nuclei passing through the crystal to the lowest 
order of perturbation theory. The higher order effects we consider in Sec. 7 where the cross 
section of the process is calculated. We finally conclude with a brief summary in Sec. 8. 
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2 Coherent elastic scattering and absorption. 



The interatomic distances in crystal, a, are large, compared to the Thomas-Fermi screening 

o 1 / o o 

radius r , a ~ 3 — 5A ^> r = r\BZ 2 ~ O.lA, where Z 2 is the atomic number of the 
target atom and r B is the Bohr radius [3]. The relevant impact parameters, b, satisfy the 
condition b <C a and the amplitudes of scattering by different atomic chains parallel to a 
given crystallographic axis are incoherent. 

The amplitude of small- angle scattering of charged particle (charge Z{) by a linear chain 
of N identical atoms in the eikonal approximation reads [1] 



Ffi(q) = l £j rf 2 bexp(.qb)(M/ / ({r j })|l-5(b,s 1 ,...s A )|^({r j })), (1) 

where \J/j and \I/ f are the initial and final state wave functions of the crystal and q is the 
2D-vector of the momentum transfer. The incident particle momentum p is assumed to be 
large enough to satisfy the condition of applicability of the straight paths approximation, 
p/q 2 ^> aN. The latter condition insures the coherence of interactions with different atoms. 

The elastic scattering corresponds to % — f and the brackets ( ) signify that an average 
is to be taken over all configurations of atoms in the ground state, 



(*({r J })|l-5(b,s 1 ,...s N )|*({r J })> 



J rf 3 ri ...rf 3 r N |^({rj})| 



N 



exp(i^x(//|b- Sj |)) 



(2) 



In (2) the total scattering phase is the sum of the phase shifts contributed by the individual 
atoms. The positions of the N atoms which make up the target are defined by the 3D- 
vectors rj, j = 1, N. The 2D- vectors Sj are the projections of these vectors on the impact 
parameter plane. We neglect all position correlations of the atoms and describe the ground 
state of crystal by the wave function |\&) such that 

l*({rj})| 2 = nWui)| 2 , (3) 
j'=i 



where the 3-D-vectors u,- are defined by = (j — l)a + u^, j = 1, TV, a = (0, 0, a) and 

u i = ( s i' z i)- 

From eq. (2) it follows that 



Fu(q) = F(q) = ip j bdbJ (qb) {l - (exp^/ifo)]^} . (4) 

Hereafter, Jq,i(x) and K 0jl (x) are the Bessel functions and the screened Coulomb phase shift 
function is 

X (jib) = -PK (jib), (5) 

with (3 = 2aZ\Z 2 and fx = r^ 1 . After integration over longitudinal variables {zj} followed 
by the azimuthal integration the term (exp(i%)) takes the form 
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(exp(ix)) = J d 2 sp(s) exp[?x(/i|b - s|)] 

= exp(— VL 2 b 2 ) \ dxexp(-x) 
Jo 

xI (2bQ^)exp[-i(3Ko{fi^/Q)]. (6) 

The 2D-vector s, describes the position of the target atom in the impact parameter plane. 
The one-particle probability distribution p(s) is as follows 

p(s) = J dz\^(s, z)\ 2 = (Q 2 /n) exp (-fiV) . (7) 

For the most commonly studied elements at room temperature the ratio p/Q varies in a 
wide range, from p/Q ~ 0.1 to p/Q ~ 1 [3]. Consider first the region of small impact 
parameters including b < 1/Vl. 2 For b < 1/2Q only small s, such that lis < 1, contribute. 
One can put then K (ps) ~ log(l//xs) and integrate over s, 

(expfe))^ (Q lP e W (-n 2 b 2 ) 

x / dxx lP/2 &v(-x)I {Wly/x) 
Jo 



, iP\*( iP.,. o2,2 



In eq.(8) 



^ . , a z a(a + 1) z 2 . . 

*( a , M= l + -- + ^- + ... (9) 



is the confluent hyper-geometric function and <3>(a; 6; z) = exp(z)<&(& — a; 6; —2). 
From eq. (8) it follows that 



|(exp(ix))| 6=0 

where the identity 



Tip 



2 sinh(7r^/2) 



1/2 

(10) 



has been used. In the weak coupling regime, (5 <C 1, 

1 



exp(ix)>| 6=0 ^l-^((x 2 >-(x> 2 ) 



and 



2 Notice that the smallness of the ratio r A /u 2 ~ 10~ 5 — 1CP 6 , where ta is the nuclear radius and u = 
is the amplitude of lattice thermal vibrations, allows one to neglect the nuclear interactions of the projectile 
up to N <~ 10 5 . As we shall see, the absorption effect which we are interested in enters the game at much 
smaller N. 
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(x 2 ) - {x? 



24 



(12) 



For > 1 



l(exp(ix))| 6=0 ~ ^7r/?exp(-7r/?/4). (13) 

Therefore, at small impact parameters, b < 1/2Q, the intensity of outgoing nuclear waves as 
a function of N exhibits the exponential attenuation. In terms of the unitarity cuts of the 
elastic scattering amplitude the imaginary part of the phase shift function comes from the 
cuts through the multi-photon projectile-atom blocks as shown in Fig. la. The account of 
diagrams like that of Fig. lb which allows cuts only between projectile-atom blocks, gives 
the scattering matrix of the form exp(iiV(x)) and affects only the overall real phase of the 
amplitude. 




a) 



b) 



Figure 1: Example of the relevant multiple scattering diagrams to order /3 3 . The unitarity 
cut of the elastic amplitude la) which contributes to the absorption is shown by crosses. The 
diagram lb) allows cuts only between the projectile- atom blocks and does not contribute to 
the absorption effect. 

The absorption effect becomes weaker toward the region of large impact parameters 

b>i/2n, 



|(exp(ix)>r ^ |(exp(ix))ir=o 



N 



i + + 



(14) 



For still larger b, b 3> 1/2Q, making use of the asymptotic form Iq(z) ~ (2nz) 1 / 2 exp(z) 
yields 



(exp(ix)) ^ 2^ 



sds 
y/nbs 



exp[-fi 2 (6 - s) 2 ] exp[ix(fJ>s)]. 



(15) 



To evaluate the integral (15) expand x(/ us ) i n powers of (s — b), 

X{ns) ~ x(vb)+u;(s-b). 
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If the frequency u, 

is small compared to fl, 



u<^n, (16) 
(x 2 ) - (x) 2 = 



cu 2 



2fi 2 



and 



(exp(?x)) ~ exp [ix — iojb] 

o \fnbs 



xfi [ S f S exp[— il 2 (b — s) 2 } exp[ius 
Jo \/Tibs 



~ exp(ix) exp[-o//41r]. (17) 

The condition (16) is satisfied if b ^> (3/Q. For the impact parameters from the region 
/3/Q <^ b <^ l//iwe may write u ~ (3/b and for larger b such that b ^> /i -1 , 



w -^y^ ex p(-^)- ( 18 ) 

From the consideration presented above it follows that the absorption effect in the elastic 
scattering amplitude is especially strong for impact parameters 

For b <^b a the atomic chain acts like an opaque "black" disc. Certainly the value of this 
finding differs for different observables and for different processes proceeding at different 
impact parameters. The only thing which is worth noticing here is the representation of the 
scattering matrix in the form 

(S(b))~e W ^iN x -^Pj-, (20) 
The eq.(20) supplemented with the observation that 

(5(6))^(7r^/ 2 exp(-^; b (21) 
simplifies all further calculations greatly. 



6 



3 The elastic cross section. 



Integrating once by parts reduces F(q) to the form convenient for evaluation of the total 
cross section, 



= roo M6Ji( 6)(i ^ (i )>( (i )>(JV -i)_ (22) 

q Jo 



ip/iN 

At small impact parameters, b <C 1/2//, 



xrf^iUf^l;^. (23) 



Because of multiple scatterings only large impact parameters b may contribute to F(q) at 
large iV and small q. Hence, 



F(q) ~ 6d6Ji(g6) [ix' - uu'/2n 2 

x exp(iiVx) exp(-A^cu 2 /4fi 2 ). (24) 



where the explicite form of (exp(i%)) at large 6, eq.(17), has been used. For large b, u 2 oc 
exp(— 2 fib) and diminishes with growing b much faster than the phase shift function 
which is oc exp(— fib). One can see that the leading contribution to the elastic scattering 
amplitude (24) comes from 

b ~ fi' 1 ^ > b a , 

where 

f = \og(PN) 

and for large iV the second term in square brackets in (24) is small compared to the first 
one. Then, for 

and £ >> 1, the steepest descent from the saddle-point 



bo = »~ 1 [£ + wr/2] 



(25) 



in eq.(24) yields 



F(q) ~ -^-Mqbo). (26) 
9 

The effect of lattice thermal vibrations at small g appears to be marginal and reduces to the 
factor exp(/x 2 / 4f2 2 N) in (26) which is irrelevant at large N. The amplitude F(q) at q — > 
coincides with the elastic scattering amplitude given by the SL approximation [2]. 

If <? ^ <?o the stationary phase approximation gives the elastic scattering amplitude of 
the form 
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where 77 = log(/i(3N/q) ^> 1. The account of the lattice thermal vibrations insures the 
convergence of the integral for the coherent elastic scattering cross section, 

which for £ ^> 1 is simply 



ffd « ^e 2 - (29) 
A* 

4 The quasi-elastic cross section. 

In this paper we focus on the coherent nucleus-atom interactions. The incoherent process of 
ionization of the target atom is suppressed by the factor ~ Z^ 1 . Then, the inelastic process 
which via unitarity gives rise to the attenuation of elastic amplitude is the process of the 
quasi-elastic scattering (Fig. la). Its cross section is the quantity [1] 

P 2 ^ = EI^*(9)I 2 -I^(?)l 2 (3°) 



d 2 q 



f 



where the sum extends over all final states of crystal, in which no particle production takes 
place. The closure relation then yields 

^ = i/^ 2 b'exp[iq(b-bO] 

-(eM^b)}) N (^PHx*^b')]) N } (31) 

and 

a Qel = Jd 2 b{l-\(ex P Mfib)])\ 2N }. (32) 

In the SL approximation | (exp[i% (//&)]) | = 1 and aq e[ = 0. From (20) and the discussion of 
the absorption radius, b a , presented above it follows that for £ 3> 1 

l-|(exp[*x(/i&)]>r *9(2b a -b) (33) 

and 

a QA » 7v(2b a ) 2 » ^ 2 (34) 
/i 
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5 The total cross section. 



From eq.(26) by means of the optical theorem we find the total cross section 

47T 27T 

a tot = — JmF(O) « — £ 2 . (35) 
p p 2 

Thus, we conclude that at high energy and in the large- N limit 

Oel ~ OQel ~ T^tot- (36) 

6 Coulomb excitation of ultra-relativistic particles and 
nuclei in crystal channel. The excitation cross sec- 
tion to the lowest order. The Born approximation. 

Consider now the process of the coherent Coulomb excitation of ultra-relativistic particles 
and nuclei passing through the crystal. This way of the experimental study of rare processes 
has been proposed in [4, 5, 6, 7, 8, 9, 10]. 

The ultra-relativistic projectile-nucleus (the mass number A, the charge Z\ and the four- 
momentum p) moving along a crystal axis undergoes a correlated series of soft collisions 
which give rise to diagonal (A — > A, A* — > A*) and off-diagonal (A — > A*, A* — > A) 
transitions. 

In [9, 4, 5] it has been proposed to study the electric dipole transition in 19 F, the 
excitation of the state \J W = 1/2 ) from the ground state |l/2 + ). The phenomenological 
matrix element of the transition l/2 + — > l/2~ is [14] 

M = -du{p') lb (qe - eq) u(p), (37) 

where both u(p') and u(p) are bispinors of initial and final states of the projectile, d is 
the transition dipole moment and e is the photon polarization vector. The transverse and 
longitudinal components the 4-vector p — p' are denoted by q and k, respectively. In what 
follows q = |q|.The only phenomenological parameter of the problem is the dipole moment 
d. The measured life-time of the 110 KeV level 19 F(l/2~) is r = (0.853 ± 0.010) x 10~ 9 
sec [13] and the dipole moment of the l/2 + — > 1/2" transition, determined from the width 
of the level 19 F(l/2~) is d ~ 5 x 10~ 8 KeV -1 [14]. Then, first, because of large value of r 
the decay of excited state inside the target crystal can be safely neglected and, second, due 
to the smallness of d, the excitation amplitude is much smaller than the elastic Coulomb 
amplitude for all q up to q ~ \fAnaZi/d and can be considered as a perturbation. Thus, 
the multi-channel problem reduces to the one-channel one. 

The high-energy helicity-flip Born amplitude of the transition l/2 + — > l/2~ in collision 
of the projectile-nucleus with iV bound atoms in crystal reads 

where a = (<Ji,a2,cr 3 ) is the Pauli spin vector, {ai,(jj} = 2Sij and the amplitude we are 
constructed is to be regarded as an operator which transforms the initial helicity state of 
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the projectile into its final state. In the denominator of eq.(38) A 2 = ^ 2 + k 2 . In the Glauber 
approximation the longitudinal momentum transfer which determines the coherency length, 



l c ~ k 1 , reads [11] 



MAE 



K = 



P 



(39) 



where M is the mass of projectile and AE is the excitation energy 3 . 
The structure factor of crystal, S(k), to the first order in g is 



S(k) = exp 



4fi 2 



sin(/tiVa/2) 
sin(/ta/2) 



(40) 



If the projectile momentum satisfies the resonance condition [4, 5, 9, 7] 

MAE 2nn 



n 



P 



0, 1, 2. 



(41) 



S(k) ~ N. Then, to the first order in g and to the zero order in (3 (Born approximation) 
the cross section of the coherent excitation of the projectile in scattering on a chain of N 
atoms in crystal is 



g 2 N 2 

47T 



log 1 + 




2Q 2 



X 2 + 29? 



(42) 



where g = y/47iadZ 2 . The central idea of [4, 5, 6, 7, 9, 10] based on the Born approximation 
is that the transition rate can be enhanced substantially due to coherency of interactions 
which is assumed to sustain over the large distance scale. The law a ex oc N 2 is expected 
to hold true up to the crystal thicknesses N = L/a ~ 10 5 — 10 6 in tungsten target. In 
[10] the Born approximation for the coherent excitation of S + in high-energy proton-crystal 
interactions p~y — > S + has been assumed to be valid up to iV ~ 10 8 . However, the account 
of the initial and final state Coulomb interactions dramatically changes the dependence of 
a ex on N. For instance, at iV = 2 the excitation amplitude is of the following form 



i4 2) (q) = - / d 2 bexp(*qb)(^exp(* X )>(exp(* X )> 

71 J 



(43) 



The first one of two bracketed factors in eq.(43) corresponds to the nuclear excitation am- 
plitude in scattering on a bound atom. It differs from the excitation amplitude of the Born 
approximation, /^(b), by the multiplicative phase factor which is due to the initial and 
final state multiple Coulomb scattering. At small impact parameters, b < l/2f2, 



(/,?exp(i x )> ~ 5(k) 



exp 



2 %2 



(44) 



For large b 



3 The Fresnel corrections to the eikonal approximation which are neglected here become important at 
large N or at large q diminish the coherency length and bring about an additional suppression of coherent 
processes [12]. 
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(f e B x exp(i X )> = S(k)£- J d 2 sp(s)(cT(n h - n s )) 
xAXi(A|n b - n s |) exp [ix (n\n b - n s |)] 

~ S^j-^n^XK^Xb) exp(* X ) exp , (45) 

where n b = b/|b|, n s = s/ jsj and b > As far as for small 6 

l(^exp(^))| 2 ocfi 2 6 2 
and for large impact parameters, b > 1//J,, 

bK 2 (fib) oc exp(— 2/i6) 



the cross section 



ff S = /d 2 b|</*expfo)}n<expfo)}| 



= ^'(1-^J <«) 

is dominated by 6 ~ 1/2//. For the diamond crystal fi/Q ~ 0.16 [3]. Hence, u 2 /2Q 2 ~ 
2f3 2 n 2 /n 2 ~ 1/20. This estimate shows that even for the diamond crystal target the Born 
approximation is irrelevant already at iV > 10. 

7 The multiple scattering effects and absorption in the 
coherent Coulomb excitation processes 

The transition amplitude on a chain of iV identical atoms including all the multi-photon 
t-channel exchanges reads 

FeM) = - I rf 2 bexp(^qb)(^expfe))(expfe)) Ar - 1 (47) 

71 J 

Because of both the multiple scattering effect and absorption only large impact parameters, 
b ^> may contribute to F ex (q). Then,the evaluation of F ex (q) reads 



2 

2 



FeM » F^(«)(^n q ) f°° bdbMqb) 

xXK^Xb) exp(iiVx) exp(-Nco 2 /An 2 ), (48) 



where n q = q/|q|. The contribution of the domain q < q = to the excitation cross 
section can be neglected as far as F ex oc q in this region. If q ^> q and £ ^> 1, the stationary 
phase approximation gives the coherent excitation amplitude of the form 
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We see that the helicity-flip dynamics removes the factor 1/q from the elastic amplitude 
(27) thus making the UV-regularization of the excitation cross section indispensable. The 
latter is evaluated as, 



7T r 2 

&ex = J d( l \ F < 
g 2 N l-8 

8tt 



ex 



a. 



|2 




(50) 



where C = 7 A A 2 r(A), 7 = 2ft 2 //3V, A = A//x and 5 = A - 1 ~ /t 2 /2 / u 2 < 1. In (50) we 
put simply S(n) = N. Thus, the account of multiple scatterings and absorption turns the 
Born approximation cross section a ex oc iV 2 into a ex oc N 1 " 6 logiV . In the limit of p — > 00 
and £ — > 0, 

^-^Tlog^-J (51) 

The dependence of o" e:E on iV differs from that of the fully unitarized elastic cross section, 
a e i oc log 2 N. The reason is that in a ex we sum the eikonal diagrams to all orders in (3 but 
only to the first order in g. Such a procedure of unitarization is, of course, incomplete, but 
this is of no importance for practical purposes since the smallness of d 2 Q 2 makes the next 
to leading order terms negligibly small up to iV ~ aZ 2 /5Vl 2 d 2 ~ 10 12 . 



8 Summary 

The main goal we pursued in this paper is a consistent description of the coherent Coulomb 
excitation of ultra-relativistic particles and nuclei passing through the aligned crystal. We 
started with the discussion of the elastic scattering and found that the account of the 
lattice thermal vibrations within the Glauber multiple scattering theory gives rise to the 
strong absorption effect. The radius of the absorption region in the impact parameter space 
appeared to grow logarithmically with growing crystal thickness. We derive convenient 
representation for the scattering matrix with absorption and calculate the coherent elastic 
and the incoherent quasi-elastic cross sections. The suppression of scattering amplitudes in 
the absorption region is shown to serve as a natural UV regulator and enables consistent 
calculation of the cross section of the coherent nuclear excitation, a ex . The dependence of 
a ex on the crystal thickness is found. The multiple scattering effects are shown to become 
numerically important already at iV > 1 thus leaving no room for the Born approximation 
widely used in early analyses of the problem. 

Acknowledgments: Thanks are due to N.N. Nikolaev for useful comments. 
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